The shortest time and/or the shortest path strategies 
in a CA FF pedestrian dynamics model 
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Abstract — This paper deals with a mathematical model of 
a pedestrian movement. A stochastic cellular automata (CA) 
approach is used here. The Floor Field (FF) model is a basis 
model. FF models imply that virtual people follow the shortest 
path strategy. But people are followed by a strategy of the 
shortest time as well. This paper is focused on how to math- 
ematically formalize and implement to a model these features 
of the pedestrian movement. Some results of a simulation are 
presented. 

Keywords -Cellular automata; pedestrian dynamics; transi- 
tion probabilities; artificial intelligence 

I. Introduction 

A stochastic cellular automata (CA) model of pedestrian 
flow is considered here. Our model takes inspiration from 
stochastic floor field (FF) CA model |2| that provides 
pedestrians with a map which "shows" the shortest distance 
from current position to a target. In this paper we focus on 
mathematical formalizing and implementation to the model 
such behavioral aspects of decision making process as: while 
moving people follow at least two strategies — the shortest 
path and the shortest time. Strategies may vary, cooperate 
and compete depending on current position. 

This is a next attempt to extend basis FF model towards 
a behavioral aspect making more flexible/realistic decision 
making process and improve simulation of individual and 
collective dynamics of people flow. 

II. Statement of the problem 

The space (plane) is known and sampled into cells 
40cm X 40cm which can either be empty or occupied by 
one pedestrian (particle) only [2|. Cells may be occupied by 
walls and other nonmovable obstacles. So space is presented 
by 2 matrixes: 

1, cell is occupied by a pedestrian; 
0, cell is empty. 



1, cell (i, j) is occupied by an obstacle; 
0, cell is empty. 



A Static Floor Field (SFF) S is used in the model. Field 
S coincides with the sampled space. A value of each Sij 
saves the shortest distance from cell to a nearest exit; 
i.e., S increases radially from exit cells where Sij are zero. 
It doesn't evolve with time and isn't changed by the presence 
of the particles. One can consider S" as a map that pedestrians 
use to move to the nearest exit. 

Starting people positions are known. A target point for 
each pedestrian is the nearest exit. Each particle can move 
to one of four its next-neighbor cells or to stay in present 
cell (the von Neumann neighborhood) at each discrete time 

step fig. m i.e., V,nax = 1. 
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Figure 1. Target cells for a pedestrian in the next time step (2)- 

Generally speaking, a direction for each particle at each 
time step is random and determined in accordance with a 
transition probabilities distribution (and transition rules). 

So a main problem is to determine "right" transition 
probabilities (and transition rules). 

III. Solution 

A. Update rules 

A typical scheme for stochastic CA models is used here. 
There is step of some preliminary calculations. Then at each 
time step transition probabilities are calculated, and direction 
is chosen. If there are more then one candidates to one 
cell a conflict resolution procedure is applied, and then a 
simultaneous transition of all particles is made. 



In our case the preliminary step includes calculations of 
SFF S. Each cell Sij saves shortest discreet distance to 
the nearest exit. The unit of such distance is a number 
of steps. To calculate the field S (and only here) we 
admit diagonal transitions and consider that a vertical and 
horizontal movement to the nearest cell has a length of 1; a 
length of a diagonal movement to the nearest cell is ^/2. (It's 
clear that movement through a corner of walls or coUums is 
forbidden and around movement is admitted in such cases 
only.) It is made a discreet distance more close to continuous 
one. 

Probabilities to move from cell to each of four the 
nearest cells are calculated in the following way: 

_ Pi-lJ _ Pi,j + 1 

* '■^ Norrrii j ' ^'"'^ Norrrii j ' 

Pi+l.J = ^TT ,Pl,j-l = T7 , (1) 

Norrriij iM orrriij 

where Norm^j ^ pi_^ j + Pi+i,j +Pi,j-i- 

Moreover 

Pi-i.j = 0, Pi,j+i = 0, Pi+i,j = 0, Pi,j-i = (2) 
only if 

Wi-ij = 1, Wij + i = 1, Wi+i,j = 1, = 1 

(3) 

correspondingly. 

A probability to stay at present cell isn't calculated 
directly. But decision rules are organized in a way that 
such opportunity may be realized, and a people patience 
is reproduced by this means. 

Decisions rules are the following [SI: 

1) If Norrriij — then motion is forbidden, otherwise 
a target cell {l,m)* is chosen randomly using the 
transition probabilities. 

2) a) If Norrriij 7^ and (1 — /;*,„) = 1 then a target 

cell {l,m)* is fixed, 
b) If Norm.^j ^ and (1 - /*„) = then the cell 
{l,m)* is not available for moving and a "people 
patience" can be realized. To do it probabilities 
of the cell {l,m)* and all other occupied the 
nearest neighbors are given to an opportunity 
not to leave the present position. A target cell is 
randomly chosen again among empty neighbors 
and the present position. 

3) Whenever two or more pedestrians have the same 
target cell, the movement of all involved pedestrians 
is denied with the probability fi, i.e. all pedestrians 
remain at their places f2l. One of the candidates moves 
to the desired cell with the probability 1 — /i . A 
pedestrian that is allowed to move is chosen randomly. 

4) Pedestrians that are allowed to move perform their 
motion to the target cell. 



5) Pedestrians that stand in exit cells are removed from 
the room. 

These rules are applied to all particles at the same time; 
i.e., parallel update is used. 

B. How to calculate probability? 

Mostly in this paper we focus on transition probabilities. 
In normal situations people choose their route carefully 
(see 1 1 1 and reference therein). Pedestrians keep a certain 
distance from other people and obstacles. The more hurried 
a pedestrian is and the more tight crowd is the more smaller 
this distance is. While moving people follow at least two 
strategies — the shortest path and the shortest time. 

In FF models people move to the nearest exit, and their 
wish to move there doesn't depend on a current distance to 
the exit. From the probability view point this means that for 
each particle among all the nearest neighbor cells a neighbor 
with the smallest S should have the largest probability. So 
a main driving force for each pedestrian is to minimize SFF 
S at each time step. But in this case only a strategy of the 
shortest path is mainly realized, and a slight regard to an 
avoidance of congestions is supposed. This is not reaUstic 
for people movement. 

An idea to improve a dynamics in a FF model is to 
introduce an environment analyzer in a probability formula. 
It should decrease an influence of a short path strategy and 
increase the possibility to move to a direction with favorable 
conditions for a moving. This will provide some kind of 
"trade off" between two main strategies. 

In this paper we introduce a revised idea of the environ- 
ment analyzer Q and make an attempt to mathematically 
formalize a complex decision making process that people 
do choosing their path — while moving their strategies 
may vary: cooperate, coincide and compete depending on 
a current position and an environment; i.e., depending on a 
place and time. 

At first let us present a probability formula and later 
we are discussing it in details. For example,the transition 
probability to move from a cell (i, j) to the up neighbor is: 

n- , — A^^^ APeople ^wall ( _ \ 

Pi-l.J ~ ^i~lj^i-lj ^i-lj\^ Wi-l^j). (^) 

Here 

• ^fj'ij — exp {kgASi^ij) — the main driven force: 

1) '^Si—ij — Sij Si— ij', 

2) fcs > — a sensitivity parameter (model param- 
eter) that can be interpreted as the knowledge of 
the shortest way to the destination point, or as 
a wish to move to the destination point, ks = 
means that pedestrians don't use information from 
the SFF S and move randomly. The higher ks is 
the more directed is movement of pedestrians. 

As far as, SFF depict direct distance from each cell to 
the nearest exit then ASi-ij > if cell (i — 1, j) is 



closer to exit than current the cell ASi-ij < 

if the current cell is closer And ASt-ij = if cells 

and (z — l,j) are equidistant to the exit. 
In contrast to other authors that deal with the FF model 
(e.g., f2], f4l, fS"], f?!) and use pure values of the field 
S in the probability formula we propose to use only 
AS'i-i.j. From a mathematical view point it is the same 
but computationally this trick has a great advantage. 
Values of SFF may be too high (it depends on a size of 
the space), and exp {ksSi^ij) is uncomputable. This 
is a significant restriction of that models. At the same 



time < AS*, 



< 1, and problem of computing 



^i-i^j is absent; 
A^lf^ = exp(-fcpA-i,,(r*_i,,)) — a factor that 
takes into account a people density in the direction: 



1) r* 



a distance to a the nearest obstacle in 



this direction (r*_^ ^ < r); 

2) r > — a "visibility" radius (a model parameter) 
is a maximal distance (number of cells) at which 
the pedestrian can look through to collect infor- 
mation about the density and possible obstacles 
(but not pedestrians); 

3) density < A-i,j(r*_ij) < 1, if all r*_i ^. cells 
are empty in this direction then Di^i j{r*_^ j) ~ 
0, if all r*_^ J cells are occupied by people in this 
direction then Di-ij{r*_i j) = 1. We estimate 
density by using idea of the kernel Rosenblat- 
Parzen's ^density estimate, and 



A-i.,(r^i,) 



m — 1 ^ '-^ ^ 



i—mj 



were 



$(z) 



(0.335 - 0.067(z)2) 4.4742, \z\ < V5; 

0; \z\ > V5, 

(5) 

4) kp > ks — a people sensitivity parameter (a 
model parameter) determines an influence of the 
people density. The higher fcp is the more pro- 
nounced is the strategy of the shortest path. 

Awall 

= exp - %^)i(A5',_ij - max AS-.j)) - 

a factor that takes into account walls and obstacles: 

1) k'iv > ks — a wall sensitivity parameter (a model 
parameter) determines an influence of walls and 
obstacles; 

2) max AS'i.j — 

= max{AS'i_ij, ASij+i, ASi+i^j, ASi^j-i}, 
'o, 0<O, 
1 otherwise. 



An idea of the function i(A5,;_ij — max AS'^ j) 
goes from a fact that people avoid obstacles only 
moving towards a destination point. But if people 
take detours (that means not minimizing the SFF) 
approaching to obstacles is not avoiding. 
< NOTE that only walls and obstacles turn the probability 
to "zero". 

Probabilities to move from cell to each of four 

neighbors are: 

Pi-i,j = exp[fcsAS'i„ij - fcpA-ijK*_ij)- 



' ^(1 - ^^)l(A5,_i,j - max AS,,j)J (1 - 



(6) 



Pij+i = exp[A;5AS'jj+i - /cpAj+i (''ij+i)- 
W(l ""-^ — )l(AS'ij+i - max AS'ij)] (1 - Wjj+i); 



(7) 



Vi+\,j = exp[A;sAS'j+ij - fcpA+i j (j-r+ij)- 
-kw{\ ^ — ^)i(AS'j+ij - max AS'jj)] (1 -Wi+ij); 



(8) 



Vi,j-\ = exp[ksASi^j^i - /cpAj-i«j-i)- 



-kwi^- 



)i(A5ij_i - max AS'ij)] (1 -mij_i); 



(9) 



In dSll-® a product ^p^p'^^A'"'^" is the environment 
analyzer that deals with people and walls. Parameters kp 
and kw allow to tune sensitivity of the model to the people 
density and the approaching to obstacles correspondingly. 
And as far as < AS* < 1, < D{r*) < 1 and 
< 1 — ^ < 1 both parameters shouldn't be less then 
ks- The term A^"-^'- is only to avoid obstacles ahead, we 
will not discuss it here and let kw = ks- 

The following the shortest path strategy means to take 
detours around high density regions if it is possible. The 
term Ap^°p^^ works as a reduction of the main driving force 
(that provides the shortest path strategy), and probability of 
detours becomes higher. The higher kp > ks is the more 
pronounced the shortest time strategy is. Note that the low 
people density makes influence of AP'^"^^^ small, and the 
probability of the shortest path strategy increases for the 
particle. 

IV. Simulations 

Here we present some simulation results to demonstrate 
that our idea works. We use one space and compare 2 sets 
of parameters. Size of space is 14.8to x 13.2™ (37 cells 
X 33 cells) with one exit (2.0m). Recall that the space is 
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Figure 3. Evacuation for 300 people, fcg = kw = 4, r = 10, kp = 6. 



sampled into cells of size 40cm x 40cm which can either be 
empty or occupied by one pedestrian only. The static field 
S is presented in fig. 2a. In fig. 2b are stating positions of 
particles. They move towards exit with v — v„iax — 1- 

Here we don't present some quantity results and only 
demonstrate a quality difference of the flow dynamics for 
2 sets of model parameters for the model presented. 

The first set of parameters is ks = kw = 4, fcp = 6, 
r = 10. The second set is ks = kw = A, kp = 18, r = 10. 
A following moving condition are reproduced by both sets - 
pedestrians know a way to the exit very well, they want go 
to the exit (it is determined by ks), a visibility is good (r), 
attitude to walls is "loyal" (kw — ks). Only one parameter 
varies here, it's fcp. 

In the first case (fcp = 6) a prevailing moving strategy is 
the shortest path. Fig. 3 presents an evacuation in different 
moments for this case. 

The other set of parameters ks = kw — 4:, kp — 18, 
r = 10 (see fig. 4) allows to realize both strategies depending 
on conditions. Recall that the term ^P<=^°p'f= in (|6])-(|9]) works 
only if the people density D{r*) > 0, and it reduces 



Figure 4. Evacuation for 300 people, ks = k^y = 4, r = 10, kp = 18. 



the probability of the shortest path strategy depending on 
density's value. 

V. Conclusion 

Figures 3-4 show a great difference in the flow dynamics 
that obtained by following only one movement strategy and 
by "keeping in mind" both strategies at a time. The case of 
fcp = 18, i.e., when both strategies of the shortest path and 
the shortest time are well pronounced, gives a more realistic 
shape of flow. A model dynamics proper needs a careful 
investigation and it is go on. A necessity of the fcp spatial 
adaptation is already clear. 
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